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I.  Introduction 


In  a  recent  publication,1  a  straightforward,  numerically  stable  method 
for  performing  electro-optic  effect  calculations  was  presented.  Simple 
analytic  procedures  were  developed  for  calculating  the  principal  dielectric 
axes  and  refractive  indices  of  an  electro-optic  crystal  of  any  crystal  class 
subject  to  an  external  electric  field  applied  in  a  general  direction.  Simple 
formulas  were  also  developed  for  obtaining  the  two  allowed  eigen- 

polarizations  and  their  corresponding  refractive  indices  for  a  general 
direction  of  phase  propagation.  However,  in  addition  to  an  external 
electric  field,  the  optical  properties  of  a  crystal  may  be  affected  by  other 
influences,  internal  or  external,  such  as  natural  optical  activity,  an 
internal  or  applied  magnetic  field,  stress,  and  others.  The  optical 

properties  and  the  induced  changes  in  them  may  be  described  by  the  relative 
permittivity  tensor  [ c ]  or  its  inverse,  the  impermeability  tensor  [*7],  where 
[ ]  =  [cf1  *  [l/nz]  and  n  is  the  index  of  refraction. 

The  linear  and  quadratic  electro-optic  effects  induce  changes  in  the 
linear  birefringence  of  a  crystal.  These  effects  may  be  represented  as 
symmetric  perturbations  to  the  impermeability  tensor.  The  allowed 
polarizations  remain  linear,  regardless  of  the  direction  of  propagation. 
However,  some  electro-optic  crystals  such  as  bismuth  silicon  oxide  are  also 
optically  active,  thereby  exhibiting  natural  reciprocal  circular 
birefringence.  This  physical  effect  is  manifested  as  a  rotation  of  the 

linear  polarization  of  the  light  upon  transmission  through  the  crystal. 
Correspondingly,  the  eigen-polarizations  are  no  longer  linear  but  are,  in 
general,  elliptical.  Rotation  of  linearly  polarized  light  may  also  be 
induced  by  an  external  field.  For  example,  an  applied  magnetic  field  may 
induce  nonreciprocal  circular  birefringence  (Faraday  rotation) . 


Furthermore,  an  external  electric  field  may  induce  reciprocal  circular 
birefringence  (electrogyration  effect).  Media  exhibiting  circular 
birefringence,  whether  natural  or  induced,  are  referred  to  as  gyrotropic . 
The  mechanisms  giving  rise  to  gyrotropy  may  be  represented  as  imaginary 
antisymmetric  (Hermitian)  perturbations  to  the  impermeability  tensor. 

The  two  general  questions  to  be  addressed  here  are  as  follows:  1)  Given 
a  crystal  that  is  linear  birefringent  (natural  or  induced)  and/or  gyrotropic 
(natural  or  induced) ,  what  are  the  principal  refractive  indices  and 

principal  dielectric  axes  of  the  crystal?  2)  What  are  the  two  eigenstates 
(i.e.,  phase  velocity  indices  and  corresponding  eigen-polarizations)  for  a 
given  direction  of  light  propagation?  The  other  given  conditions  are  that 
the  crystal  is  lossless,  linear,  and  homogeneous.  The  answers  to  the  above 
questions  can  be  obtained  by  first  determining  the  eigenvalues  and 

eigenvectors  of  the  perturbed  impermeability  tensor.  In  addition,  since  the 
eigen-polarizations  are  now  elliptical  in  general,  three  pieces  of 

information  are  required  to  describe  each  state  of  polarization:  azimuth 

angle,  ellipticity,  and  handedness. 

Similar  problems  with  the  electromagnetic  description  of  gyrotropy  have 
been  addressed  in  the  literature  in  terms  of  macroscopic  theory  and 
quantum  mechanics.9  11  As  an  example  of  the  former,  a  method  to  calculate 
the  eigenstates  of  a  naturally  optically  active,  electro-optic  crystal  by 
diagonalizing  the  coupled-wave  equations  was  presented  by  Yariv  and 

Lotspeich.3  As  an  example  of  the  latter,  a  method  using  quantum 

electrodynamics  to  determine  the  eigenstates  in  optically  active,  linear 
birefringent  crystals  was  presented  by  Eimerl.10  The  approach  presented  in 
this  paper  is  based  on  the  macroscopic  properties  of  the  crystal.  The 
procedure  introduced  here  employs  an  extension  of  the  general  Jacobi  method, 
a  known,  very  accurate,  stable,  and  simple  numerical  routine.  Also,  to 


obtain  a  full  description  of  the  eigen-polarizations,  a  complex  polarization 
variable  is  used  that  maps  the  eigenvectors  of  the  transverse  impermeability 
tensor  into  a  complex  polarization  plane.  The  advantages  of  this  method 
over  the  others  are  that  (1)  no  assumptions  are  required;  (2)  it  applies  to 
any  crystal  class;  (3)  it  applies  to  any  field  direction;  (4)  it  applies  to 
any  direction  of  light  propagation;  and  (5)  it  is  numerically  stable, 
accurate,  and  straightforward. 

The  constitutive  equation  will  first  be  used  to  describe  the  optical 
properties  of  a  crystal.  To  provide  a  geometric  interpretation  of  linear 
and  circular  birefringence,  the  index  ellipsoid  and  gyration  surface  will  be 
reviewed.  Next,  the  procedure  to  address  the  given  problem  will  be 
introduced,  which  includes  a  brief  overview  of  the  eigenvalue/eigenvector 
problem  for  Hercitian  matrices,  followed  by  a  description  of  the  extended 
Jacobi  method  introduced  in  the  present  work  and  of  the  complex  polarization 
variable  used  to  obtain  the  eigenstates  for  a  given  direction  of  phase 
propagation.  Finally,  bismuth  silicon  oxide,  an  optically  active,  electro¬ 
optic,  electrogyratory  material  is  analyzed  to  illustrate  the  simplicity  and 
accuracy  of  the  method.  Throughout  this  paper  boldface  will  be  used  to 
denote  a  vector  and  [•]  to  denote  a  matrix. 

II.  Constitutive  Equation 

The  material  properties  (principal  permittivities  or  refractive 
indices)  of  the  crystal  are  represented  by  the  constitutive  relation  D  - 
[c]E,  where  [«]  is  the  permittivity  tensor  of  the  medium.  Disturbances  to 
the  optical  properties,  which  are  typically  very  small  in  magnitude,  may  be 
described  by  this  tensor.  They  are  commonly  expressed  in  terms  of  the 
inverse  of  the  permittivity  tensor,  [ e  ]  ~ 1  -  1/c  [c]"1  -  l/co[l/nz],  where 
«o  is  the  permittivity  of  free  space. 


For  a  homogeneous,  lossless,  and  nongyrotropic  medium  the  permittivity 
(impermeability)  tensor  has  only  real  components.  Moreover,  it  is  symmetric 
for  all  crystal  classes  and  for  any  selection  of  the  dielectric  axes.12' 14 
Therefore,  it  may  be  diagonalized,  and  in  principal  coordinates,  the 
constitutive  equation  is 
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with  the  principal  permittivities  on  the  diagonal.  The  symmetry  of  f<] 
guarantees  that  a  diagonal  form  exists  given  a  correct  choice  of  three 
perpendicular  principal  axes  (x,y,z)  with  respect  to  the  crystallographic 
axes . 

If  the  medium  is  gyrotropic,  the  constitutive  equation  may  be  written 
as13'15 

D  -  [ e ]E  +  ieoG  x  E  -  [e]'E,  (2) 

where  [«]  is  the  symmetric  unperturbed  permittivity  tensor,  i  is  J-l ,  and 
G  is  the  gyration  vector  which  is  uniquely  defined  for  the  mechanism 
producing  the  circular  birefringence.  The  vector  cross  product  G  x  E  in 
Eq.  (2)  may  be  represented  as  the  product  of  an  antisymmetric  tensor  [G] 
with  the  vector  E.  Thus, 


D  -  {[«]  +  ieo [G] }E  -  [c]'E.  (3) 

Therefore,  the  permittivity  tensor  [c]'  is  now  clearly  Hermitian  as  indeed 
it  must  be  due  to  thermodynamic  arguments.13’14 

The  constitutive  equation  may  also  be  written  as 


(A) 


E  -  [«]'-xD  -  M'D  -  (l/f0)(  [»?]  -  i[»f]  [GJ  [«,)  ID. 


From  Eq.  (4),  the  antisymmetric  (imaginary)  part  of  [17] '  is 
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An  important  point  is  that  the  imaginary  part  of  [1?]'  has  no  effect  on  the 
principal  dielectric  axes  and  indices  of  the  crystal.  The  real  part  of  [17]' 
is  real  and  symmetric. 

A  number  of  the  various  types  of  influences,  both  natural  and  induced, 
on  the  optical  properties  are  now  described. 

A.  Dielectric  Properties  with  No  Fields  Applied 

1.  Natural  Linear  Birefringence 

The  optical  symmetry  of  a  crystal  is  represented  by  the  permittivity 
(impermeability)  tensor.12'14  If  all  diagonal  elements  are  equal,  then  the 
crystal  is  isotropic.  If  e*  -  ey  *  tz  ,  then  the  crystal  is  uniaxial.  If  * 
ty  r*  tz  *  ,  the  crystal  is  biaxial.  Therefore,  both  uniaxial  and  biaxial 

crystals  exhibit  natural  linear  birefringence;  e.g.,  ez  -  ey  **  0. 

2.  Natural  Optical  Activity 

In  general,  the  macroscopic  properties  of  a  medium  depend  on  the 
temporal  and/or  spatial  variation  of  the  electromagnetic  field.  For  the 

case  of  natural  optical  activity,  the  properties  are  influenced  by  spatial 
dispersion,  the  dependence  of  [e]  on  the  magnitude  and  direction  of  k  at 
fixed  frequency.4,13'19  The  macroscopic  dipole  moment  per  unit  volume  of  the 
medium  at  a  given  point  depends  on  the  field  at  and  near  that  point.  In  the 
optical  frequency  range,  the  effects  of  spatial  dispersion,  in  general,  are 
small  and  are  characterized  by  the  first  power  of  a/ A  («1),  where  a  is  on 
the  order  of  the  lattice  constant  and  A  is  the  wavelength  of  the  light  in 


the  medium.  Therefore,  to  first  order,  the  constitutive  equation  is1*'15.17'18 

Di  -  £ijEj  +  -  £'ijEi  (6) 

and 

£'u  “  *  £ij(w)  +  kJf  (7) 

where  tiJ(w)  is  the  permittivity  tensor  without  optical  activity  and  7iJ^(u>) 
is  a  third-rank  real  antisymmetric  tensor  in  the  indices  i  and  j ,  resulting 
in  «*  being  Hermitian.  The  second  term  can  also  be  represented  by  G  x  E, 
where  G  is  the  gyration  (axial)  vector.13’15 

The  displacement  vector  D  rotates  in  a  helical  fashion  about  the 
wavevector  k,  so  G  is  parallel  to  k  and  the  components  of  G  are 
functions  of  the  direction  cosines  of  k.  The  sense  of  rotation  bears  a 
fixed  relation  to  the  direction  of  propagation  as  shown  in  Fig.  la.  If  a 
linearly  polarized  input  light  wave  is  transmitted  through  an  optically 
active  crystal  and  then  is  reflected  back  through  the  crystal,  the  net 
rotation  of  the  polarization  is  zero.  Therefore,  natural  optical  activity 
is  a  reciprocal  effect.*'13 

By  defining  the  direction  of  k  in  spherical  coordinates  (4>k  ,  0k  )  ,  the 
components  and  magnitude  of  G  are13,17 


G 
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-  Gsin0k  cos<f>k ,  G  -  Gsinffk  sin^k ,  G  -  Gcos#k 

7  *  (8) 

G  -  glxsin20k  cos2^k  +  g22  sin2  \  sin2  <f>k  +  g33cosz0k 

+  2g12sin20k  sin^k  cos^k  +  2g13sin0k  cos$k  cos<&k 

+  2g23sin0k  cos &k  sin^k  , 


where  gAJ  are  the  components  of  the  gyration  tensor.  Only  noncentrosymmetric 
crystals  can  have  natural  optical  activity.  Table  1  provides  a  summary  of 
the  gyration  tensors  for  all  crystal  classes  that  exhibit  optical  activity.13 


B.  Dielectric  Properties  with  External  Fields  Applied 
1.  Linear  Electro-Optic  Effect 

An  electric  field  applied  in  an  arbitrary  direction  to  a  crystal 
lacking  a  center  of  symmetry  produces  a  change  in  the  coefficients  (l/n2)i  due 
to  the  linear  electro-optic  effect  according  to 

A(l/n2)i  -  Sr  E  i  -  1 . 6  (9) 

J  j  -  x ,y , z  -  1,2,3 

where  riJ  is  the  ijth  element  of  the  linear  electro-optic  tensor  in  reduced- 
subscript  notation.13,20  In  matrix  form  Eq.  (9)  is 


A(l/n2)' 
A(l/n2)2 
A(l/n2  ) 
A(l/n2)< 
A(l/n2 )5 
A(l/n2 )6 


(10) 


For  a  noncentrosymmetric  crystal,  the  new  impermeability  tensor  in  the 
presence  of  an  applied  electric  field  is,  in  general,  no  longer  diagonal  in 
the  original  axes  system.  The  perturbed  impermeability  tensor  is 


'1/n,2  +  A(l/n2)1  A(l/n2)6  A(l/n2)5 

[ 1/n2 ] *  -  A(l/n2 )g  l/n22  +  A(l/n2)2  A(l/n2)A  (11) 

A(l/n2 )5  A(l/n2 )A  l/n32  +  A(l/n2 ),.  . 


However,  the  field- induced  perturbations  are  real  and  symmetric,  so  the 
symmetry  of  the  tensor  is  not  disturbed.  These  changes  to  [1/n2]  have  the 
effect  of  changing  the  principal  axes  and  indices  of  the  crystal.  The 
electro-optic  tensor  for  all  crystal  classes  is  summarized  in  numerous 
texts.13,15  All  optically  active  crystals  are  also  electro-optic,  but  the 


converse  is  not  true. 


2.  Electrogyration  Effect 

An  applied  electric  field  may  not  only  induce  linear  birefringence,  but 
in  many  cases,21  if  may  also  induce  circular  birefringence  through  the 
electrogyration  (EG)  effect.  Considering  only  first-order  spatial 
dispersion,  the  gyration  tensor  is21'23 

+  ^ijk^k  ’ 

where  gAJ  is  the  gyration  tensor  with  no  applied  field  and  piJk  is  the  third- 
rank  electrogyration  tensor,  which  has  the  same  symmetry  as  the  electro¬ 
optic  tensor.  The  second- rank  tensor  alters  gAJ  in  the  same  manner  as 

Eq.  (9)  alters  the  impermeability  tensor  in  Eq.  (11).  The  net  result  is  a 
change  in  the  magnitude  of  the  gyration  vector  G  as  found  by  Eq.  (8)  and 
therefore,  a  change  in  the  specific  rotation  (polarization  rotation  per  unit 
thickness)  of  the  medium.  The  EG  effect  may  be  viewed  as  induced  optical 
activity,  a  reciprocal  effect.  This  perturbation  changes  the  impermeability 
tensor  in  the  same  manner  as  optical  activity,  that  is,  producing  Hermitian 
off-diagonal  elements  as  in  Eq.  (4)  rather  than  the  real  symmetric  off- 
diagonal  elements  of  Eq.  (11). 

3.  Faraday  Rotation 

In  a  Faraday  active  medium,  the  macroscopic  properties  are  influenced 
by  frequency  dispersion  caused  by  an  applied  magnetic  f ield . * ’ 8 • 1 3 ~ 1  * ' 2  * " 2 7 
The  field  creates  a  relative  shift  between  the  phase  velocity  indices  of 
refraction  of  the  two  allowed  eigen-polarizations,  inducing  circular 
birefringence.  The  constitutive  equation  is 

D  -  [  c  ] E  +  ifoV>B  x  E,  (13) 


where  [f]  is  the  unperturbed  permittivity  tensor  and  r/>  is  the 


magnetogyration  constant  of  the  medium.  The  real  part  of  the  impermeability 
tensor  is  symmetric  in  B,  and  the  imaginary  part  is  antisymmetric  in  B.14 
Since  only  first-order  effects  are  considered,  the  real  part  of  the  tensor 
remains  unchanged.  In  this  case,  the  imaginary  part  is  represented  by  the 
gyration  vector  that  is  proportional  to  B;  i.e.,  G  -  ^B.  Thus,  the 

components  of  G  are  functions  of  the  direction  cosines  of  B  (tB,0B),  rather 
than  those  of  k.  That  is, 

Gx  -  ^Bsin?B  cosdB  ,  Gy  -  VBsin0B  sin^B  ,  Gz  -  V>Bcos0E  .  (I4) 

The  magnitude  of  G  is  merely  the  product  of  ^  and  the  magnitude  of  B.  For 
Faraday  rotation,  the  sense  of  rotation  bears  a  fixed  relation  to  B,  as 
shown  in  Fig.  lb.  The  eigen-polarizations  are  preserved  upon  reflection  so 
that  the  net  rotation  is  doubled.  Thus,  Faraday  rotation  is  a  nonreciprocal 
effect.4 ■ 11 • 13 

C.  Combined  Effects 

If,  for  example,  a  lossless,  biaxial,  and  optically  active  crystal  is 
subject  to  an  applied  electric  field,  then  the  impermeability  tensor  is 
altered  by  Hermitian  perturbations  and  is5 ,J3~ 14,24 


[1/n2]'  -  [r,]' 


(15) 


where  Vii,  i  **  j  ,  are  complex  and  r?ii  are  real.  As  stated  previously,  the 
imaginary  parts  of  the  off-diagonal  elements  do  not  affect  the  principal 
axes  or  indices  of  the  crystal.  They  affect  only  the  state  of  the  allowed 
polarizations  and  the  phase  velocity  indices.  The  eigen-polarizations  are 
now  elliptical,  in  general,  rather  than  linear  as  with  the  electro-optic 
effect.  An  Hermitian  matrix  may  be  represented  by  a  quadratic  surface  in 


complex  space.  In  real  Cartesian  space,  however,  only  the  real  part  of 
fij]',  which  is  symmetric,  contributes  to  a  quadratic  surface  (ellipsoid). 
Geometric  surfaces  which  represent  optical  properties  of  the  crystal  are 
discussed  in  the  next  section. 

III.  Geometric  Approach 
A.  Index  Ellipsoid 

The  index  ellipsoid  is  a  construct  whose  geometric  characteristics 
represent  the  phase  velocities  and  the  directions  of  electric  displacement 
vibration  of  the  two  allowed  plane  waves  corresponding  to  a  given  optical 
wavevector  direction  in  a  crystal.  The  general  index  ellipsoid  for  an 
optically  biaxial  crystal  is  expressed  in  Cartesian  coordinates  as12"13 

(x2/nx2)  +  (y2/ny2)  +  (z2/n,2)  -  1  (16) 

where  nx ,  ny  ,  and  nz  are  the  principal  refractive  indices  of  the  crystal. 
Since  the  permittivity  (Impermeability)  tensor  is  positive  definite,  the 
surface  is  always  an  ellipsoid.  If  nx  -  ny ,  the  surface  becomes  an  ellipsoid 
of  revolution,  representing  a  uniaxial  crystal.  An  isotropic  crystal  is 
represented  by  a  sphere  (degenerate  ellipsoid)  with  the  principal  axes 
having  equal  length.  These  surfaces  are  shown  in  Fig.  1  of  Ref.  1.  Also, 
shown  are  the  optic  axes  for  each  crystal  symmetry. 

The  eigenstates  for  an  arbitrary  direction  of  propagation  in  a  crystal 
are  found  in  the  elliptical  cross-section  perpendicular  to  k  which  passes 
through  the  origin  of  the  index  ellipsoid,  as  shown  in  Fig.  2.  If  the 
optical  properties  are  not  disturbed,  the  major  and  minor  axes  of  the  cross- 
section  ellipse  represent  the  two  allowed  orthogonal  linear  vibration 
directions  of  D  (eigen-polarizations)  for  that  particular  direction  of 
propagation.  Tie  lengths  of  these  axes  correspond  to  the  the  respective 


phase  velocity  indices  of  the  allowed  polarizations.  They  are,  therefore, 
referred  to  as  the  "fast"  and  "slow"  axes.  As  with  the  principal  axes  and 
indices  of  the  crystal,  these  eigenstates,  in  general,  are  affected  by  real 
symmetric  perturbations,  (for  example,  the  electro-optic  effect),  to  the 
impermeability  tensor.  However,  the  antisymmetric  perturbations  affect  only 
the  eigenstates  (the  eigen-polarizations  and  phase  velocity  indices) .  The 
major  axes  of  the  cross-section  ellipse  correspond  to  the  major  axes  of  the 
allowed  polarizations  in  this  case. 

B.  Gyration  Surface 

The  geometry  of  the  index  ellipsoid  provides  only  partial  information 
about  the  eigenstates  if  the  allowed  polarizations  are  not  linear.  In  this 
case  the  ellipsoid  can  be  used  to  determine  only  the  orientation  of  the 
allowed  elliptical  polarizations  but  nothing  on  the  properties  of  optical 
rotation  for  a  given  wavevector  direction.  The  gyration  surface,  however, 
may  be  used  to  illustrate  the  directional  dependence  of  optical  rotation  in 
a  gyrotropic  crystal.  For  optically  active  crystals  the  surface  is 

constructed  from  the  gyration  tensor  [g]  in  the  same  way  that  the  index 
ellipsoid  is  constructed  from  the  impermeability  tensor.17  That  is,  the 
distance  from  the  origin  to  any  point  on  the  surface  is  given  by 

G  -  gij  \  kj ,  (17) 

where  k±  and  k^  are  the  direction  cosines  of  the  wavevector  k.  Equation  (17) 
is  the  directional  magnitude  of  the  gyration  vector,  G  -  | G f ,  as  given  by 
Eq.  (8).  Since  the  gyration  tensor  is  not  necessarily  positive  definite, 
the  surface  may  have  a  variety  of  forms.  Shubnikov17  provides  a  complete  set 
of  all  possible  surfaces  for  isotropic,  uniaxial,  and  biaxial  crystal 
classes  which  exhibit  optical  activity.  For  example,  Fig.  3  is  the  gyration 


surface  for  right-handed  quartz  (class  32)  which  is  positive  uniaxial  and 
optically  active.  The  first  two  diagonal  elements,  gx  -  g22 ,  are  negative, 
and  g33  is  positive.  The  surface  is  given  by  G  -  -|g1 x |sin2 0ko  +  g33cos20ko, 
where  0ko  is  the  angle  between  the  optic  axis  and  k.  Therefore,  optical 
rotation  along  the  optic  axis  is  right-handed  and  is  denoted  by  the  white 
surface.  However,  propagation  perpendicular  to  the  optic  axis  gives  optical 
rotation  in  the  opposite  direction  and  is  denoted  by  the  dark  surface. 
There  is  no  optical  rotation  (G  -  0)  for  propagation  directed  -56°  from  the 
optic  axis,  as  determined  by  the  measured  quantity  g11/g33  ~  0.45,  and  the 
eigen-polarizations  are  linear. 

For  Faraday  active  crystals  a  gyration  surface  may  also  be  constructed. 

As  stated  before,  the  sense  of  optical  rotation  bears  a  fixed  relation  to 
the  direction  of  the  applied  magnetic  field  B.  Maximum  rotation  is 
achieved  for  propagation  parallel  and  antiparallel  to  the  applied  magnetic 
field.  If  propagation  is  in  a  direction  inclined  to  B,  the  degree  of 
optical  rotation  will  decrease  as  cos#kB,  where  0kB  is  the  angle  between  k 
and  B.141  Therefore,  Fig.  4  is  a  representation  of  the  gyration  surface  for 
this  type  of  crystal. 

IV.  Analysis 

The  problem  is  to  determine  the  principal  axes  and  indices  of  the 
crystal  and  the  two  allowed  orthogonal  eigenstates  of  polarization  (DJl  and 
D2 )  and  the  corresponding  phase  velocity  indices  (n1  and  n2 )  for  a  general 
wavevector  direction  and  a  general  external  field.  The  solution  to  the 

problem  lies  in  determining  the  eigenvalues  and  eigenvectors  of  Eq .  (15). 

The  method  chosen  to  address  this  problem  involves  diagonalizing  the  matrix 
using  an  extension  of  the  general  Jacobi  method.  For  Hermitian  matrices, 
the  eigenvalues  are  real.  The  eigenvectors,  on  the  other  hand,  are  complex 


in  general.  Thus,  additional  information  is  required  to  describe  the 
general  state  of  the  eigen-polarizations  (the  complex  eigenvectors).  For 
linear  polarization,  only  the  orientation  (azimuth  angle)  in  the  plane 
transverse  to  k  is  required.  For  elliptical  polarization,  the  ellipticity 
and  handedness  as  well  as  the  azimuth  angle  must  be  determined. 
Furthermore,  for  linear  polarization,  D  oscillates  in  a  plane  in  one  fixed 
direction  perpendicular  to  k.  If  D  is  elliptically  polarized,  then  it  no 
longer  oscillates  in  a  plane  but  rather  propagates  in  a  flattened  helix 
about  k,  as  shown  in  Fig.  5.  A  full  description  of  the  eigen-polarizations 
is  obtained  by  using  a  complex  polarization  variable. 

A.  Principal  Axes  and  Principal  Refractive  Indices 

A  detailed  description  of  a  procedure  to  determine  the  new  principal 
axes  and  indices  of  a  crystal  subject  to  real  symmetric  perturbations  is 
given  in  Ref.  1.  That  approach  employs  the  general  Jacobi  method,  and  it  is 
known  to  be  an  accurate,  numerically  stable,  and  simple  routine  for 
diagonalizing  real  symmetric  matrices.  The  routine  is  iterative,  and  it 

involves  the  calculation  of  an  elementary  plane  rotation  angle  at  each  step 
to  zero  the  largest  off-diagonal  element.  Also,  the  method  allows  for 
consistent  labeling  of  the  new  axes;  a  global  rotation  axis  and  a  global 
rotation  angle  are  calculated  from  the  resulting  cumulative  orthogonal 
transformation  matrix.  To  find  the  principal  axes  and  indices  of  a  crystal 
with  Hermitian  perturbations  to  its  optical  properties,  as  given  by  Eq . 

(15),  the  same  procedure  is  performed  but  only  on  the  real  part  of  the 

matrix.  The  authors  refer  the  reader  to  Ref.  1  for  this  part  of  the 

analysis . 

B.  Eigenstates  of  Polarization  and  Phase  Velocity  Indices 

The  eigenvalue  problem  for  Hermitian  matrices  is  addressed  by  a  unitary 


transformation,  [a][»j3'[a]B  -  [A],  where  [a]  is  the  unitary  transformation 
matrix  ( [  a  ] H  -  [a]"1),  [rj]'  is  an  Hermitian  matrix,  [A]  is  the  resulting 
diagonal  matrix  of  real  eigenvalues,  and  H  denotes  complex  conjugate 
transpose.  As  suggested  by  Wilkinson,28  a  form  of  the  unitary  matrix  that 
can  be  used  is 

la]  —  T  cos4>  eiBsin3>l 

[-e~ 1 8  s  in4>  cos$  J.  (18) 

This  matrix  has  the  effect  of  transforming  a  system  from  Cartesian 
coordinates  to  a  complex  (helical)  coordinate  system.  There  are  two 
parameters,  $  and  B,  to  determine.  For  real  symmetric  matrices  B  -  0,  and  $> 
represents  the  elementary  Cartesian  plane  rotation  angle.  Using  the 
expression  for  a  unitary  transformation,  a  set  of  relationships  was  derived 
{with  the  unitary  matrix  of  Eq.  (18)),  which  results  in  a  version  of  the 
general  Jacobi  method  extended  to  Hermitian  matrices.  The  two  unknowns  $>  and 
B  are  calculated  at  each  iteration  step.  With  additional  algebra,  the 
simple  expressions  in  the  Appendix  were  obtained  for  updating  the  elements 
of  [rj]'  as  the  diagonalization  process  proceeds.  The  parameter  B  was 
determined  to  be  the  argument  of  the  off-diagonal  element  i  **  j  ,  and  is 

the  required  value  for  zeroing  that  element  with  a  rotation  in  the  (i,j) 
complex  plane.  These  formulas  reduce  to  those  of  the  general  Jacobi  method 
for  real  symmetric  matrices.1  The  Jacobi  method  using  these  modified 
formulas  was  programmed  and  tested.  The  results  for  several  test  matrices 
were  found  with  virtually  zero  error.  These  results  were  often  more 

accurate  than  those  obtained  with  the  commercial  IMSL30  routine,  EIGCH. 

To  find  the  eigenstates  for  a  specific  direction  of  propagation  k, 
a  real  orthogonal  transformation  must  first  be  performed  on  the  Hermitian 
matrix  [»/]'  to  place  the  problem  in  a  coordinate  system  of  k.  The 


wavevector  direction  is  specified  by  the  spherical  coordinates  (^k  ,  0k  )  in  the 


original  (x,y,z)  coordinate  system.  A  new  coordinate  system  (x' ' ,y' ' ,z* ' ) 
is  defined  with  z''  parallel  to  k  and  x''  lying  in  the  (z,z'')  plane.  The 
transformation  to  the  (x' ' ,y' ' ,z' ' )  system  is  produced  first  by  a  rotation 
^k  about  the  z-axis  followed  by  a  rotation  6 k  about  y' '  as  shown  in  Fig.  6 
This  transformation  is  described  by 


x'  '  -  x  cos#k  cos <£k  +  y  cos#k  sin^k  -  z  sin#k 

y' '  -  -x  sin£k  +  y  cos<f>k  (19) 

z''  -  x  sin#k  cos^k  +  y  sin#k  sin^k  +  z  cos 8k  . 

The  transformed  tensor  [r/]''  in  the  (x'’,y'')  plane  is 


[•?]" 


(20) 


where 

»7XX"  "  ^xx  cos2^k  +  ’Jyy  sin2^k)  cosz*k  +  f7zisin2^k  + 

2 cos28k  cos^k  sin^k  -  2cos0ksin0k  (i)jzrcos <£k  +  »7yir  sin^k), 

»?  ’  '  -  n  sin2^.  +  tj  cos2i  -  2f?  cos0.  sin^.  ,  (21) 

»?xy"  -  (»?yy  -  ,?JCX)cosek  cos^k  sin^k  +  cos*k  (»!„  cos2^k  -  f;*xy  sin2^k ) 

+  sin0k  (fj*xz  sin^k  -  >7*yz  cos^k )  -  >7*yx  "  . 

and  r7iJr  denotes  the  real  part  of  »?1J  .  Therefore,  the  Hermitian  property  is 
preserved.  The  third  row  and  column  can  be  neglected,  since  the  vibration 
direction  of  the  eigen-polarizations  is  transverse  to  k.  Next,  the  phase 
velocity  indices  of  refraction  are  easily  determined  using  the  relationships 
in  the  Appendix  (with  a  slight  modification)  to  diagonalize  the  2x2  matrix 
of  Eq.  (20).  The  required  rotation  in  the  complex  (x'',y'*)  plane  is 


*  -  bTan* 1  ( 2  sgnO?^  "  >  b„  *  '  |/0»,x  '  '  -  1 


xy 


-yy 


)) 


(22) 


where  |-|  denotes  magnitude  and  sgn  (.flxyt'')  is  the  sign  of  rlxyi' '  •  The  phase 
velocity  indices  are 


(^^'cos2}  +  t;yy''sin24  +  2  sgn(r?xy  r  '  '  )  |  ??xy  '  '  j  cos$  sin#)"15 
{fjix"sin24>  +  *7yy  '  'cos24>  -  2  sgn(fjxy  r  "  )  |  nxy  "  |  cos4>  sin?-}'15 


(23) 


The  corresponding  states  of  polarization  are  just  the  rows  of  the  2x2 
unitary  transformation  matrix  of  Eq.  (18): 


Dx  -  [cos4  eiBsin$]B 
D2  -  [-e*iBsin'I>  cos4>]B  , 


(24) 


where  B  -  Arg (»;'').  These  states  are  the  left  eigenvectors  of  [rj]'';  the 
right  eigenvectors  are  the  complex  conjugate  transpose  of  Eq.  (24)  and  are 


D:  -  j  cos® 


e~  iBsin$>J  , 


D2  -  |-e  sin$ 


The  orthogonality  relation,  Dj-D2*  -  0,  is  satisfied.  From  Eq .  (23)  if  nx  < 
n2 ,  as  shown  in  Fig.  7,  then  D:  is  the  "fast"  wave  and  D2  is  the  "slow"  wave. 
If  n:  >  n2 ,  then  the  "fast"  and  "slow"  waves  are  reversed. 

Now  the  full  description  of  these  eigen-polarizations  must  be 
determined;  i.e.,  azimuth  angle,  ellipticity,  and  handedness.  They  have  the 
same  ellipticity  but  opposite  senses  of  rotation  with  orthogonal  major  axes. 
First,  the  azimuth  angle  or  direction  of  the  major  axes  for  and  D2  is 
easily  obtained  from  Eq.  (22)  as  follows.  The  orthogonal  transformation 
performed  on  [>;]'  can  be  visualized  in  terms  of  the  index  ellipsoid.  The 
cross-section  ellipse  transverse  to  k  is  found  by  taking  only  the  real  part 
of  the  transformed  tensor  [»?]''.  giving  X'  ' T Re {  [ rj ]  '  '  )X'  '  -  1 ,  or 


'x' 


r)  '  'y' 

'  77  J 


,2 


2r,  "x-y  -  1, 


(26) 


where  Re  {[»?]''}  is  the  real  part  of  [  *j  ]  *  *  and  X''  -  [x*  '  y'']r.  The  azimuth 
angle  f}^  for  Dj  is  the  angle  in  the  real  Cartesian  (x",y")  plane  required 
to  diagonalize  Re([r/]''}.  It  is 


“  Wan"1  [2i»xyr' '/(•?„"  -  V")]-  (27) 

This  expression  is  the  same  as  Eq .  (22)  if  the  numerator  2sgn (*?xyr")  |  y  '  '  | 

is  replaced  by  2 »7xyr''-  The  azimuth  angle  for  D2  is  just  P2  “  +  K /^-  The 

angles  ^  and  $2  define  the  orthogonal  semi-axes  directions  (x'''  and  y' ' '  in 
Fig.  7)  of  the  cross -section  ellipse  of  the  index  ellipsoid.  If  I '*) 
were  zero,  the  lengths  of  the  semi -axes  would  correspond  to  the  phase 
velocity  indices  for  the  two  linear  polarizations.  These  indices  would  be 
calculated  from  Eq.  (23)  with  /9  rather  than  4>. 

The  ellipticity  and  handedness  of  Dj  and  D£  can  be  found  through  the 
use  of  a  complex  polarization  variable  (CPV)  *.31  These  eigen-polarizations 
are  in  the  form  of  a  two-component  Cartesian  Jones  vector  orthogonal  to  k: 
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(28) 


The  form  of  the  CPV  is  then  x  “  rexp(iAS),  where  r  -  | Dy /Dx |  and  A 5  - 

6^-8*.  It  performs  a  bilinear  transformation  from  the  complex  (x'',y'')  plane 
to  polarization  space,  shown  in  Fig.  8,  where  points  in  the  complex  x  plane 
represent  polarization  states.  This  plane,  in  fact,  maps  directly  onto  the 
Poincare'  sphere.  The  equator  of  the  sphere  is  the  horizontal  real  axis  of 
linear  polarizations,  and  the  north  and  south  poles  are  the  points  R  and  L 
of  the  right-  and  left-circular  polarizations,  respectively.  From  D: ,  the 


CPV  is 


e'  iBtan4>, 


(29) 


*1 

and  from  the  orthogonal  state  D2 , 

X2 - 1/Xj*  “  -e'iBcot*.  (30) 

Relationships  between  the  CPV  xx  and  the  azimuth  angle  f. 9  and  ellipticity 

angle  are  given  by  Azzam  and  Bashara3 1  as 


tan2^  -  2Re(x1 )/ (1  -  \X,  |2) 

sin2^  -  2Im(x1)/(l  +  |xj2)- 


(31) 


In  terms  of  the  impermeability  tensor  elements  and  the  rotation  angle  *,  the 
azimuth  angle  0^  was  given  in  Eq.  (27),  and  the  ellipticity  angle  £  was 
derived  to  be  sin2^1  -  -{Im(fjxy  ' ' ) / 1 »7X y  |  )  sin2*  -  -sinB  sin2*  or 

-  -h  sin*1(sinB  sin24>)  .  (32) 

The  ellipticity  is  given  by  tan£: .  Furthermore,  the  relative  amplitude  of 
the  orthogonal  components  of  D:  is  defined  as  r:  -  tan  *,  and  the  relative 
phase  is  A ^ - Arg(i7xy") - B. 

The  handedness  of  the  polarization  is  determined  by  the  sign  of  the 
relative  phase.  If  AS  is  >  0,  then  the  elliptical  polarization  is  left- 
handed.  If  AS  <  0,  then  the  polarization  is  right-handed.  And  if  AS  -  0  or 
m*,  m  -  1,2,3...,  then  the  polarization  is  linear. 

For  the  orthogonal  polarization  state,  the  parameters  £2 ,  e2,  r2 ,  and  A S2 
for  x2  ®re31 

P2  -  +  w/2,  {2  -  ,  r2  -  cot*,  A52  - -  +  AS1.  (33) 

The  orthogonality  condition  for  Xj  and  Xz  is  X x  X2*  *  -1.  Therefore,  in  the 


process  of  diagonalizing  a  Hermitian  matrix,  the  relative  phase  A 6  and 
ellipticity  angle  £  are  changing  iteratively.  In  the  case  of  a  real 
symmetric  matrix,  the  process  is  interpreted  as  a  rigid  body  rotation.1 

The  advantages  of  this  method  include  the  following:  (1)  It  is  accurate 
and  stable  for  all  crystal  classes;  (2)  Stable  orthonormal  eigenvectors  are 
found  simultaneously  with  the  eigenvalues;  (3)  The  unitary  transformation 
matrix  is  easily  determined  by  elements  of  the  perturbed  impermeability 
tensor;  (4)  All  descriptive  information  about  the  eigen-polarizations  of  a 
crystal  for  a  given  k  is  obtained  from  the  unitary  transformation  matrix 
with  simple  formulas;  and  (5)  The  combined  effects  of  real  and  Hermitian 
perturbations  to  the  impermeability  tensor  can  be  straightforwardly  handled. 

V.  Example:  The  Sillenite  Crystal  Class 

A  cubic  sillenite  crystal  of  class  23  is  examined  to  illustrate  the 
ease  and  accuracy  of  the  method  just  described.  This  crystal  class  is  being 
widely  investigated  for  applications  in  dynamic  real-time  holographic 
interferometry  and  spatial  light  modulation.23,32  The  class  includes  bismuth 
silicon  oxide  (BSO) ,  bismuth  germanium  oxide  (BGO) ,  and  bismuth  titanium 
oxide  (BTO) .  These  crystals  are  electro-optic,  optically  active, 
electrogyratory ,  piezoelectric,  and  elastooptic.  The  electro-optic  effect, 
optical  activity,  and  electrogyration  effect  combined  may  strongly  affect 
the  polarization  in  these  crystals.  BSO  is  examined  to  illustrate  these 
influences  on  the  eigen-polarizations. 

Two  principal  configurations  of  BSO,  both  transverse,  are  used  for 
volume  holography.  One  of  them  is  shown  in  Fig.  9. 3 2  An  external  electric 
bias  is  applied  in  the  [110]  direction,  and  the  direction  of  light 
propagation  is  [1  1  0]  or  (^k.#k)  “  (135°, 90°).  In  the  (x,y,z)  coordinate 
system,  the  unperturbed  impermeability  tensor  for  this  cubic  isotropic 


crystal  is 


[l/n,2  0  0 

[1/n2]  -  0  l/i^2  0  (34) 

0  0  l/i\,aJ. 

where  n0  is  the  principal  refractive  index  of  the  crystal.  The  index 
ellipsoid  is  a  sphere.  The  gyration  tensor  for  natural  optical  activity  is 
a  function  of  the  given  applied  electric  field  and  wavevector  directions  and 
is 


[g]'  - 


b1 

S"  *  1  Ey 


f*lEx 


lEy 
8i  i 


where  f  is  the  electrogyratory  coefficient.  The  gyration  surface  changes 
from  a  sphere  to  the  surface  {from  Eqs.  (8)  and  (17)), 


G'  -  glx  -  72  f41  E  sin$k  cos0k  (sin^k  +  cos^k ) ,  (36) 

where  E  -  E  -  -1/72  E,  and  E  is  the  magnitude  of  applied  field.  If  E  -  0, 
then  G'  -  G  -  g  ,  which  is  a  sphere.  The  perturbed  gyration  surface  is 
shown  in  Fig.  10.  Note  by  Eq.  (36)  that  for  the  given  field  ([1  1  0])  and 
wavevector  ([1  1  0])  directions,  the  effect  of  electrogyration  is  not 
present.  Therefore,  with  the  electro-optic  effect,  optical  activity,  and 
the  given  direction  of  propagation  k  taken  into  account,  the  tensor  becomes 

l/no2  0  l/72(-r41E  -  igll/noM 

111'  -  0  l/no2  l/72{-r41E  -  igll/no*) 

l/72{-r41E  +  igll/no4)  1/721-r^E  +  igxl/noM  1/n,2 

(37) 

where  Gz  -  -G  -  -1/72  g: x .  The  new  orientation  of  the  index  ellipsoid  is 
obtained  by  applying  the  general  Jacobi  method  as  described  in  Ref.  1  to  the 


real  part  of  Eq.  (37). 

The  crystal  is  now  biaxial , 

and  the 

principal 

indices  and  axes  are 

n  .  -  2.52996 

x'  -  [1/2  1/2  -1/72]1 

n  .  -  2.53 
y  * 

y'  -  [-1/72  1/72  0]T 

(y'll  k> 

(38) 

nr,  -  2.53004 

z'  -  [1/2  1/2  1/72] T . 

To  determine  the  eigenstates  of  the  crystal  for  the  given  k,  the  entire 
[17]'  tensor  of  Eq.  (37)  must  be  transformed  to  the  (x'  '  ,y '  '  , 2 '  '  )  coordinate 
system  by  Eqs .  (19).  The  resulting  2x2  matrix  in  the  (x'',y'')  plane  is 
given  by 


[»7 


ft  _ 


1/n 


L(-r41E  -  ig^/nj) 


(”r4 1 E  +  igxl/no*) 
1/n  2 


(39) 


and  x"  -  -2  ([0  0  1]),  y"  -  -1/72  x  -  1/72  y  ([1  1  0]),  and  z"  - 
-1/72  x  +  1/72  y  ([I  1  0]  ||  k)  .  The  cross-section  ellipse  is,  by  Eq.  (26), 

l/nQ2  (x"2  +  y"2)  -  2r41E  x"y"  -  1.  (40) 


The  azimuth  angle  is  found  by  placing  this  cross-section  in  principal 
coordinates,  i.e.,  by  diagonalizing  the  real  part  of  Eq.  (39).  This  angle 
is  -  HTan' 1 {-2r4 1E/(l/nQ2  -  l/nQ2))  -  -45°.  Therefore,  the  axes  of  the 
cross-section  ellipse  are  along  x'''  -  [1/2  1/2  -1/J2}T  and  y',F  - 

[-1/2  -1/2  -1/72]T.  Numerical  values32  for  the  various  parameters  are  nQ  - 

2.53  and  r41  -  4.41xl0’12  m/V  at  a  freespace  wavelength  of  -  0.6328  /im.  From 
these  numbers  the  lengths  of  the  principal  axes  are  nx,,,  -  1^  —  hn^r^E 
2.52996  and  ny ,  ,  .  -  +  Hn^r^E  -  2.53004  for  a  field  magnitude  of  106  V/m. 

The  phase  velocity  indices  are  found  by  diagonalizing  the  entire  matrix 
[i?]''  of  Eq.  (39).  The  required  rotation  angle  is 


*  -  hTan" 1  { 2  sgn(qxyr  )  |  l/d/V  -  l/t^2)}  -  -45*.  (41) 

and  |*?xy''|  -  ((r41E)2  +  (ga  j  A^*  )2  •  The  constant  g21  is  calculated  from  the 

measured  specific  rotation  of  p  -  21.4c/mm  (Ref.  32)  which  is  equal  to 

*Sn/AlV  Therefore.  gn/n0*  “  Ap/™03  "  2  6618  x  10*  -  Eiving  Nxy  "  |  “ 

2.6621x10**.  Optical  activity  dominates  the  magnitude  of  the  off-diagonal 
element  .  The  phase  velocity  indices  are  found  from  Eqs.  (23)  to  be. 


nx  -  (l/no2  +  1  *7x y  ' '  1  } “  2.52785  (fast  wave) 
n2  -  {l/nQ2  -  |»?xy  "()'**  -  2.53216.  (slow  wave) 


(42) 


The  circular  (elliptical)  birefringence  is  then  Anc  -  n2  -  n:  -  0.0043111. 

Without  the  electric  field  applied,  the  indices  are  nL ,  -  2.52785  (-  nx )  and 
n2 ,  -  2.532158  (<  n2 ) ,  giving  a  circular  birefringence  of  An. ,  -  n2 ,  -  nx ,  “ 
0.0043080.  Therefore,  the  electro-optic  effect  only  slightly  enhances 
circular  birefringence.  The  corresponding  eigenvectors  are  found  from  Eq. 
(25)  to  be 


1/72 

-(0.01656  +  iO. 99986)l//2. 


"(0.01656  -  iO . 99986) l/v 

1/72 


and  B  -  Arg(qxy'')  -  -89.0508°. 

The  CPV  is  -exp(-iB).  From  Eq.  (31),  the  azimuth  angle  is  -  -45°, 
which  agrees  with  the  angle  obtained  before.  Also  from  Eq.  (31),  the 
ellipticity  angle  is  sin2^  -  Im(x  )  “  sinB,  which  is  equal  to 

— Im(»jxy  ’  '  )/ j *7xy  '  '  |  sin2$  -  -sin B  sin24>  -  sinB  from  Eq.  (32)  for  $  -  -45°. 
Therefore,  Eqs.  (31)  and  (32)  are  consistent.  The  ellipticity  angle  is 
then  —44.52541°  and  the  ellipticity  of  the  polarization  is  tan£x  -  -0.98357° 
(almost  circular  polarization).  The  relative  amplitude  is  rt  -  1,  and  the 
relative  phase  between  Dxl  and  Dy x  is  L61  -  +89.05°.  The  corresponding  point 


on  the  complex  x  plane  in  Fig.  8  is  on  the  unit  circle  at  the  phase  angle  of 
+89.05°.  Since  ASj^  >  0,  then  Dx  is  left-handed  polarization,  and  nx  -  is 
the  corresponding  index.  Therefore,  D2  corresponds  to  the  fast  wave. 

For  the  orthogonal  polarization  J>2  ,  the  CPV  is  x2  “  exp(-iB).  The 
azimuth  angle  is  +45°,  and  the  ellipticity  angle  is  — £  or  +44.52541°  for  an 
ellipticity  of  0.98357.  The  relative  amplitude  is  r2  -  1,  and  the  relative 
phase  between  Di2  and  Dy2  is  -90.95°.  The  corresponding  point  on  the  complex  x 
plane  in  Fig.  8  is  on  the  unit  circle  at  the  phase  angle  of  -90.95°.  Since 
the  M2  <0,  D2  is  right-handed  polarization,  and  n2  -  .  Therefore,  T>2 

corresponds  to  the  slow  wave. 

Finally,  an  additional  note  is  that  the  phase  retardation,  T,  between 
Dx  and  D2  for  a  given  crystal  thickness  d  may  be  calculated  from  the 
circular  (elliptical)  birefringence  An^  .  It  is  given  by  r  -  2tt/A  Ancd  and  in 
general,  includes  the  effects  of  both  the  natural  optical  activity  and  the 
external  electric  field  (e.g.,  electro-optic  and  electrogyration  effects). 

VI.  Conclusion 

A  straightforward  systematic  procedure  for  performing  electro-optic 
effect  calculations  was  developed  and  presented  in  Ref.  1.  That  approach, 
which  employs  the  general  Jacobi  method,  can  be  used  to  analyze  propagation 
in  electro-optic  materials  in  any  crystal  class  for  an  arbitrary  electric 
field  direction  and  arbitrary  wavevector  direction.  The  properties  of  the 
impermeability  tensor  were  exploited  to  arrive  at  simple,  stable,  and 
accurate  expressions  for  determining  the  principal  axes  and  indices  of  a 
crystal  and  the  eigenstates  (phase  velocity  indices  and  eigen-polarizations) 
for  a  given  direction  of  propagation  k. 

In  this  paper  that  procedure  has  been  extended  to  gyrotropic  crystals. 
External  (or  internal)  influences,  such  as  optical  activity,  electrogyration 


effect,  and  Faraday  rotation  may  now  be  included,  singly  or  together, 
these  circular  birefringence  effects  are  more  complicated  and,  in  general, 
produce  elliptical  eigen-polarizations.  The  extended  method  requires  a 
unitary  transformation  from  a  Cartesian  coordinate  system  to  a  complex 
helical  coordinate  system  to  determine  the  eigenstates  for  a  given  direction 
of  propagation.  Using  a  unitary  matrix  suggested  by  Wilkinson,28  a  set  of 
formulas  has  been  derived  which  results  in  an  extended  version  of  the 
general  Jacobi  method  applicable  to  Hermitian  matrices.  These  relationships 
are  given  in  the  Appendix.  Furthermore,  a  complex  polarization  variable  was 
introduced  to  quantify  the  connection  between  the  elements  of  the  perturbed 
impermeability  tensor  and  the  eigen-polarizations.  This  new  procedure 
reduces  easily  to  the  less  complicated  case  where  the  electro-optic  effect 
alone  is  present. 

More  specifically,  two  questions  were  posed  and  answered  in  the  present 

work: 

1)  Given  a  crystal  that  is  linear  birefringent  (natural  or  induced) 
and/or  gyrotropic  (natural  or  induced) ,  what  are  the  principal 
refractive  indices  and  principal  dielectric  axes  cf  the  crystal? 

2)  What  are  the  eigenstates  for  an  arbitrary  direction  of  light 
propagation? 

The  step-by-step  procedure  introduced  in  this  paper  is  summarized  as 
follows : 

1)  The  general  Jacobi  method  as  described  in  Ref.  1  is  applied  to  the 
real  part  of  the  perturbed  impermeability  tensor  [ij]'  to  determine 
the  principal  indices  and  axes  of  the  crystal. 

2)  A  real  orthogonal  transformation  is  performed  on  [ij]*  to  place  the 
tensor  in  the  coordinate  system  of  k  (x*  *  ,y ' ' ,  z  *  ’ )  ,  giving  [»?]''. 
Using  formulas  from  the  Appendix,  the  eigenstates  for  the  given  k 


3) 


are  obtained. 


A)  From  the  eigenvectors  of  [ij]*'  a  complex  polarization  variable 
(CPV)  is  defined.  Using  the  CPV,  the  descriptive  properties  of 
the  eigen-polarizations,  i.e.,  azimuth  angle,  ellipticity  and 
handedness,  are  determined  in  terms  of  the  elements  of  [*?}''. 

To  provide  a  geometric  interpretation  of  linear  and  circular  birefringence, 
the  index  ellipsoid  and  gyration  surface  were  used. 

Finally,  the  sillenite  crystal  class  was  examined  to  illustrate  the 
ease  and  accuracy  of  the  extended  method.  Specifically,  bismuth  silicon 
oxide  (BSO)  was  analyzed  in  a  principal  configuration  to  show  the  effects  of 
its  natural  optical  activity  together  with  the  simultaneous  influences  of  an 
applied  electric  field  (through  the  electrogyration  and  electro-optic 
effects)  on  the  eigenstates  of  the  crystal  for  a  given  k.  Example 
numerical  results  were  presented. 


Appendix:  Complex  Plane  Rotations 


The  general  Jacobi  method  is  modified  here  for  diagonalizing  a  3x3 
Hermitian  matrix.  This  iterative  procedure  involves  a  unitary 
transformation  from  a  Cartesian  to  a  helical  coordinate  system.  As  a  quick 
reference  the  following  relationships  for  the  rotation  angle,  the  unitary 
transformation  matrix,  and  the  updated  impermeability  elements  are  provided 
for  rotation  in  each  of  the  three  complex  planes.  These  expressions  are 
derived  from  the  coordinate  transformation  law  for  second-rank  tensors.  The 
rotation  angle  in  the  complex  (x,y)  principal  plane  is  denoted  by  <f> ,  in  the 
complex  (x,z)  principal  plane  by  6,  and  in  the  complex  (y,z)  principal  plane 
by  \{> .  The  impermeability  tensor  is  represented  as  the  Hermitian  matrix  H. 

1.  Rotation  in  the  complex  (x,y)  plane. 

The  rotation  angle  <f>  required  to  zero  the  H1Z  element  is  found  by 

4>  -  hTan'  1  [  2  |  Hj  2  |  /  (Hj  2  -  E^)}.  (Al) 

This  angle  represents  a  counterclockwise  rotation  in  the  complex  (x,y) 
plane.  The  transformation  matrix  is 


[a] 


4>  ~ 


cos  <f> 

-e" 1 A  sim£ 

0 


e1  Astn<t 
cos^ 

0 


0 

0 

1 


(A2) 


where  exp(iA)  -  H12/|H  2  .  The  elements  of  H  are  updated  as  follows: 


H< 

H4 

H< 

H< 

H. 


* <f>\  1 

“  H11 

cos 2  <t> 

+  H22 

sin2<^  ■ 

+  2  |H:  2  |  cos4>siT\<f> 

'^2  2 

“  Hi: 

sin2^ 

+  H22 

cos2^  • 

-  2  | H2  2  |  cos^sinijS 

^3  3 

“  H33 

'^12 

-  <H2 

2  -  H 

1 1 > (Hi 

2/Ihi2 

I  )COS<f>siTuf> 

+  H12  (cos 2<f> 

'^1  3 

-  HX3 

cos  <t> 

+  H23(H12/|H 

x 2  \ ) s ln<f>  - 

H  '* 

3  1 

'^2  3 

-  ~H1 

3<H12 

•/|h12 

| )sin^ 

+  H23cos  <(> 

-  H  '*  . 

3  2 

(A3) 


-  sin2^)  - 


H 

nz  1 
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2.  Rotation  in  the  complex  (x,z)  plane. 

The  rotation  angle  6  required  to  zero  the  H13  element  is  found  by 

6  -  man'1[2|Hl3  |/(H1X  -  H,3)]  (AA) 

This  angle  represents  a  clockwise  rotation  in  the  complex  (x,z)  plane.  The 
transformation  matrix  is 

cos 6  0  eiBsin 6 

[a]e  0  10  (A5) 

-e'iBsin$  0  cosfl  _  . 

where  exp(iB)  -  H  3/|h  |.  The  elements  of  H  are  updated  as  follows: 


- 

H11cos2tf  +  H33sin20  +  2 | 3 | cosfisinfl 

W 

mm 

H22 

H*33# 

ms 

H11sin2^  +  H33cos20  -  2 |H1 3 | cosffsinff 

(A6) 

H*12’ 

- 

Hl2COS*  +  H23*  ^H13/IH13  -  H21- 

H*13' 

(h3 3  -  H11)(H13/|H13 I )cos0sin0  +  H13(cos2*  -  sin2*)  -  H31* 

0 

W 

- 

~h12* (h13/|h13 | )sin*  +  H23cos0  -  H32'*. 

3.  Rotation  in  the  complex  (y,z)  plane. 

The  rotation  angle  t/>  required  to  zero  the  H23  element  is  found  by 

*  “  HTan-1[2|H23|/(H22  -  H,,)].  (A7) 

This  angle  represents  a  counterclockwise  rotation  in  the  complex  (y,z) 
plane.  The  transformation  matrix  is 

1  0  0 

[a]^  “  0  cost/)  eiCs inTp  (A8) 

0  -e’iCsin\£  cost/> 


where  exp(iC)  -  H23/|H23|.  The  elements  of  H  are  updated  as  follows: 


*W  "»» 

Vz'  “  H22cos2sS  +  H3  3sin2^  +  2tH23|cosV>sin^ 

^33'  -  H22sinZ^  +  H33COs2^  "  2lH23|cos^sin^ 

^2'  -  H12COSV>  +  Hl3(H23*/lH23l>Sin*  “  H2l'* 

“  -H12(H23/iH23l>Si^  +  »13C0S^  “  H3l'* 

^23'  “  <H33  "  H22)(H23/lH23i)COS^SinV5  +  H23<CO®2*  ~  '  H32*' 


(A9) 


-  0. 
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Table  I.  Gyration  tensors  gi  ^  for  all  crystal  classes  exhibiting  natural 


optical  activity.  (Ref.  13) 


Triclinic 

1 


81 1 

81  2 

81  3 

82 1 

822 

823 

-83 1 

83  2 

83  3 

Biaxial 


Monoclinic 


2  (2  ||  y) 


m  (m  i  y) 


LSj  3 


S33J 


823 


Orthorhombic 
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FIGURE  CAPTIONS 


Fig.  1.  The  sense  of  optical  rotation  relative  to  the  direction  of 

propagation  k  for  (a)  natural  optical  activity  and  (b)  Faraday  rotation. 

Fig.  2.  The  index  ellipsoid  cross-section  (cross-hatched)  that  is  normal  to 
the  wavevector  k  and  passes  through  the  origin.  The  principal  axes  of  the 
cross-hatched  ellipse  represent  the  directions  of  the  allowed  linear 
polarizations  Dx  and  D2 .  Dx ,  D2 ,  and  k  form  an  orthogonal  triad. 

Fig.  3.  Gyration  surface  for  right-handed  quartz  (class  32).  The  white 

surface  depicts  right-handed  optical  rotation  with  the  maximum  rotation 
occurring  for  propagation  along  the  optic  axis.  The  dark  surface  depicts 
left-handed  rotation  with  maximum  rotation  along  a  direction  perpendicular 
to  the  optic  axis.  There  is  no  optical  rotation  for  propagation  -56°  from 
the  optic  axis. 

Fig.  4.  Gyration  surface  for  Faraday  active  crystals.  Maximum  optical 
rotation  occurs  for  p.  gation  parallel  and  antiparallel  to  B.  The  white 
surface  depicts  rt  nation  of  one  sense  while  the  dark  surface  depicts 
rotation  of  rhe  opposite  sense. 

Fig.  5.  The  flattened  helical  contour  of  an  elliptically  polarized 

propagating  wave  at  an  instant  of  time.  The  radial  vectors  to  the  contour 
represent  the  displacement  vector  D  (l  k). 

Fig.  6.  Orthogonal  transformation  of  the  (x,y,z)  dielectric  axes  to  the 

(x' ' ,y ' ' , z ' ' )  coordinate  system  of  the  wavevector  k  (z' ’  fl  k)  represented 
in  polar  coordinates  (^k ,  8^  ) . 

Fig.  7.  The  cross-section  ellipse  of  the  index  ellipsoid  in  the  (x'',y'') 
plane.  The  x' ' '  and  y'''  axes  represent  the  major  axes  of  the  eigen- 
polarizations  oriented  relative  to  x' '  and  y' ' .  The  two  eigenstates  have 
orthogonal  major  axes,  opposite  handedness,  and  the  same  ellipticity.  The 


wavevector  k  and  the  z' '  axis  are  normal  to  the  plane  of  the  figure. 

Fig.  8.  Cartesian  complex  plane  of  polarization.  Each  point  in  the  plane 
represents  a  polarization  state.  The  basis  states  are  the  horizontal  linear 
polarization  at  the  origin  and  the  vertical  linear  polarization  at  infinity. 
The  dashed  circle  represents  the  unit  circle  (unit  relative  amplitude).  The 
radial  line  represents  a  contour  of  constant  relative  phase  of  ir/A. 

Fig.  9.  A  principal  transverse  crystal  orientation  of  Bi12SiO20  (BSO) .  The 
external  electric  field  is  applied  in  the  [1  1  0]  direction,  and  the 
direction  of  propagation  is  along  [110].  The  (x,y,z)  coordinate  system 
represents  the  unperturbed  dielectric  axes.  The  new  coordinate  system 
resulting  from  the  electro-optic  effect  is  represented  by  (x',y',z').  The 
coordinate  system  of  the  wavevector  k  is  given  by  (x' '  ,y '  '  , z '  '  )  with  z' '  ||  k. 
Finally,  the  (x' '  '  ,y '  '  '  , z '  '  ' )  coordinate  system  with  z '  '  '  ||  z''  ||  k 
represents  the  principal  axis  coordinate  system  of  the  eigenstates  for  the 
given  k. 

Fig.  10.  Gyration  surface  for  BSO  in  the  (x,z)  or  (y,z)  principal  plane. 
The  dashed  circle  represents  the  surface  projection  with  no  applied  electric 
field,  and  therefore,  the  optical  rotation  is  invariant  with  the  wavevector 
k  direction.  The  solid  "heart -shape"  contour  depicts  the  the  surface  when 
a  field  is  applied  in  the  [110]  direction.  The  magnitude  of  radial 
vectors  from  the  origin  to  the  surface  is  a  measure  of  the  optical  rotation 
per  unit  length  which  depends  on  the  direction  of  k.  Note  that  optical 
rotation  is  not  affected  by  the  electric  field  for  propagation  in  the  (x,y) 


plane . 
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